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ABSTRACT 

We report the first weak-lensing detection of a large-scale filament funneling matter onto the 
core of the massive galaxy cluster MACSJ0717.5+3745. 

Our analysis is based on a mosaic of 18 multi-passband images obtained with the Ad- 
vanced Camera for Surveys aboard the Hubble Space Telescope, covering an area of ~ 1 x 20 
arcmin 2 . We use a weak-lensing pipeline developed for the COSMOS survey, modified for the 
analysis of galaxy clusters, to produce a weak-lensing catalogue. A mass map is then com- 
puted by applying a weak-gravitational-lensing multi-scale reconstruction technique designed 
to describe irregular mass distributions such as the one investigated here. We test the result- 
ing mass map by comparing the mass distribution inferred for the cluster core with the one 
derived from strong-lensing constraints and find excellent agreement. 

Our analysis detects the MACSJ07 17.5+3745 filament within the 3 sigma detection con- 
tour of the lensing mass reconstruction, and underlines the importance of filaments for the- 
oretical and numerical models of the mass distribution in the Cosmic Web. We measure the 
filament's projected length as ~ 4.5 hz\ Mpc, and its mean density as (2.92+0.66)xl0 8 /2 74 M 
kpc~ 2 . Combined with the redshift distribution of galaxies obtained after an extensive spectro- 
scopic follow-up in the area, we can rule out any projection effect resulting from the chance 
alignment on the sky of unrelated galaxy group-scale structures. Assuming plausible con- 
straints concerning the structure's geometry based on its galaxy velocity field, we construct a 
3D model of the large-scale filament. Within this framework, we derive the three-dimensional 
length of the filament to be 18 \c\ Mpc. The filament's deprojected density in terms of the 
critical density of the Universe is measured as (206 ± 46) x p cr j t , a value that lies at the very 
high end of the range predicted by numerical simulations. Finally, we study the distribution 
of stellar mass in the field of MACSJ0717. 5+3749 and, adopting a mean mass-to-light ratio 
(M*/Lk) of 0.73 + 0.22 and assuming a Chabrier Initial-Mass Function, measure a stellar 
mass fraction along the filament of (0.9 ± 0.2)%, consistent with previous measurements in 
the vicinity of massive clusters. 

Key words: cosmology: observations - gravitational lensing - large-scale structure of Uni- 



1 INTRODUCTION 



* E-mail: mathilde.jauzac@gmail.com (MJ) 



In a Universe dominated by Cold Dark Matter (CDM), such as 
the one parameterised by the ACDM concordance cosmology, hi- 
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erarchical structure formation causes massive galaxy clusters to 
form through a series of successive mergers of smaller clusters 
and groups of galaxies, as well as through continuous accretion 



tribution on very large scales jBond et al Jl 1 996 


; Yess & ShandarirJ 


1 19961: lAraeon-Calvo et al.l |2007|; lHahn et alj 


2007b predict that 



these processes of merging and accretion occur along preferred di- 
rections, i.e., high ly anisotropically. The result is the "cosmic web" 
dBond et alii 1996b . a spatially highly correlated structure of inter- 
connected filaments and vertices marked by massive galaxy clus- 
ters. Abundant observational support for this pict ure has been pro- 
vided by large-scale gal a xy redshift surveys ( e.g. jGeller & Huchral 
1 19891 : lYork et al.l |2000| : IColless et al] 1200 ll) showing voids sur- 
rounded and connected by filaments and sheets of galaxies. 

A variety of methods have been developed to detect filaments 
in surveys, among the m a "friends of friends" algorithm (FOF, 
iHuchr a & Geller 1982) combined with "Shapefin ders" statistics 
dSheth & Jainl l2003h: th e "Skeleton" algorithm dNovikov et"al] 
|2006j; Sousbie et al a two-dimensional technique developed 

by |Moodvetal.l i ll983l) : a nd the Smoothed Hessian Major Axis Fil- 
ament Finder (SHMAFF. lBond et alj|2010h . 

Although ubiquitous in large-scale galaxy surveys, filaments 
have proven hard to characterise physically, owing to their low 
density and the fact that the best observational candidates often 
turn out to be not primordial in nature but the result of recent 
cluster mergers. Specificall y, attempts to study th e warm-hot in- 
tergalactic medium (WHIM. Icct & Ostrikeifl999l) . resulting from 
the expected gravitational heating of the intergalactic medium in 
filaments, remain largely inconclusive because it is hard to ascer- 
tain for filaments near cluster whether spectral X-ray features orig- 
inate froni_Jhe_^ament or from past or ongoing clusters merg- 
ers jKaastra et alj 120061 : iRasmussenl 120071 : iGaleazzi etakl 120091 : 
IWilliams et alj|2010h . Som e detections appear robust as they have 
been repeatedly confirmed dFang etalj|2002ll2007l :l Williams et alj 
120071) but are based on just one X-ray line. An alternative observa- 
tional method is based on a search f or filamentary overdensities of 
galaxies relative to the background jPimbblet & Drinkwaterll2004l : 
libeling et ai1l2004l) . When conducted in 3D, i.e., including spec- 
troscopic galaxy redshifts, this method is well suited to detecting 
filament candidates. It does, however, not allow the determination 
of key physical properties unless it is supplemented by follow-up 
studies targeting the presumed WHIM and dark matter which are 
expected to constitute the vast majority of the mass of large-scale 
filaments. By contrast, weak gravitational lensing offers the tan- 
talising p ossibility of detect ing directly the total mass content of 
filaments <MeadetalJl201Ch . since the weak-lensing signal arises 
from luminous and dark matter alike, regardless of its dynamical 
state. 

Previous weak gravitational lensing studies of binary clusters 
found tentative evidence of filaments, but did n ot result in clear de- 
tections. One of the first efforts was made by IClowe et~al Ul998h 
who reported the detection of a filament apparently extending from 
the distant cluster RX J1716+67 (z = 0.81), using images obtained 
with the Keck 10m telescope and University of Hawaii (UH) 2.2m 
telescope. This filamentary structure would relate two distinct sub 
clusters detected on the mass and light maps. The detection was 
not confirmed though. Almost at the same time lKaiser et alj dl998h 
conducted a weak lensing study of the supercluster MS0302+17 
with the UH8K CCD camera on the Canada France Hawaii Tele- 
scope (CFHT). Their claimed detection of a filament in the field 
was, however, questioned on the grounds that the putative fila- 
ment overlapped with both a foreground structure as well as with 
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Figure 1. Area of our spectroscopic survey of MACSJ0717. 5+3745. Out- 
lined in red is the region covered by our Keck/DEIMOS masks; outlined 
in blue is the area observed with HST/ACS. Small circles correspond to 
objects for which redshifts were obtained; large filled circles mark cluster 
members. The black contours show the projected galaxy density (see Ma et 
al. 2008). 



gaps between CCD chips. Indeed, iGavazzi et al ] j2004l) showed 
the detection to have been spurious by means of a second study 
of MS0302+17 using the CFHT12K camera. A weak gravitational 
lensing analysis w ith MPG/ESO Wide Field Imager conducted by 
iGrav et al.1 J2002I) claimed the detection of a filament in the triple 
cluster A901/902. The candidate filament appeared to connect two 
of the clusters and was detected in both the galaxy distribution and 
in the weak-lensing mass map. However, this detection too was 
of low significance and coincided partly with a gap between two 
chips of the camera. As in the case of MS0302+ 17, a re-analysis of 
the A901/A902 comp lex using high-quality HST/ACS images by 
iHevmans et alj |2008j) failed to detect the filament and led the au- 
thors to conclude that the earlier detection was caused by residual 
PSF systematics and limitation of the KS93 mass reconstruction 
used in the study by IGrav et all J2002h. A further dete ction of a 
filament candidate was reported bv lDietrich etal.N2005h based on 
a weak gravitational lensing analysis of the close double cluster 
A222/A223. However, as in other similar cases, the proximity of 
the two clusters connected by the putative filament raises the pos- 
sibility of the latter being a merger remnant rather than primordial 
in nature. 

In this paper we describe the first weak gravitational analy- 
sis of the yery massive cluster MACSJ07 17.5+3745 (z = 0.55 ; 
lEdge et al.ll2003l : lEbeling et al.ll2004 120071 ; iMa et al.ll2008l|2009l) . 
Optical and X - ray a nalyses of the system JEbeling et al.1 12004| ; 
IMa et alj 120 08. 2009) find compelling evidence of a filamentary 
structure extending toward the South-East of the cluster core. Us- 
ing weak-lensing data to reconstruct the mass distribution in and 
around MACS J07 17.5+3745, we directly detect the reported fila- 
mentary structure in the field of MACSJ0717.5+3745. 

The paper is organized as follows. After an overview of 
the observational data in Section 2, we discuss the gravitational 
lensing data in hand in Section 3. The modeling of the mass 
using a multi-scale approach is described in Section 4. Results are 
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Table 1. Overview of the HST/ACS observations of MACSJ07 17.5+3745. *: Cluster core; observed through F555W, rather than F606W filter. 



F606W F814W 



R.A. (J2000) 


Dec (J2000) 


Date 


Exposure Time (s) 


Date 


Exposure Time (s) 


Programme ID 


07 17 32.93 


+37 45 05.4 


2004-04-02* 


4470* 


2004-04-02 


4560 


9722 


07 17 31.81 


+37 49 20.6 


2005-02-08 


1980 


2005-02-08 


4020 


10420 


07 17 20.38 


+37 47 07.5 


2005-01-27 


1980 


2005-01-27 


4020 


10420 


07 17 08.95 


+37 44 54.3 


2005-01-27 


1980 


2005-01-27 


4020 


10420 


07 17 43.23 


+37 47 03.1 


2005-01-30 


1980 


2005-01-30 


4020 


10420 


07 17 20.18 


+37 42 38.8 


2005-02-01 


1980 


2005-02-01 


4020 


10420 


07 17 54.26 


+37 44 49.3 


2005-01-27 


1980 


2005-01-27 


4020 


10420 


07 17 42.82 


+37 42 36.3 


2005-01-24 


1980 


2005-01-25 


4020 


10420 


07 17 31.39 


+37 40 23.3 


2005-02-01 


1980 


2005-02-01 


4020 


10420 


07 18 05.46 


+37 42 33.6 


2005-02-04 


1980 


2005-02-04 


4020 


10420 


07 17 54.02 


+37 40 20.6 


2005-02-04 


1980 


2005-02-04 


4020 


10420 


07 17 42.79 


+37 38 05.7 


2005-02-05 


1980 


2005-02-05 


4020 


10420 


07 18 16.65 


+37 40 17.7 


2005-02-05 


1980 


2005-02-05 


4020 


10420 


07 18 05.22 


+37 38 04.9 


2005-02-05 


1980 


2005-02-05 


4020 


10420 


07 17 53.79 


+37 35 52.0 


2005-02-05 


1980 


2005-02-05 


4020 


10420 


07 18 27.84 


+37 38 01.9 


2005-02-08 


1980 


2005-02-08 


4020 


10420 


07 18 16.40 


+37 35 49.1 


2005-02-08 


1980 


2005-02-08 


4020 


10420 


07 18 04.97 


+37 33 36.2 


2005-02-09 


1980 


2005-02-09 


4020 


10420 



discussed in Section 5, and we present our conclusions in Section 6. 

All our results use the ACDM concordance cosmology with 
fl M = 0.3, £2 A = 0.7, and a Hubble constant H = 74 km s" 1 Mpc -1 , 
hence 1" corresponds to 6.065 kpc at the redshift of the cluster. 
Magnitudes are quoted in the AB system. 



2 OBSERVATIONS 

The MAssive Cluster Survey (MACS. libeling et alj|200lh was the 
first cluster survey to search exclusively for very massive clusters 
at moderate to high redshift. Covering over 20,000 deg 2 and us- 
ing dedicated optical follow-up observations to identify faint X- 
ray sources detected in the ROSAT All-Sky Survey, MACS com- 
piled a sample of over 120 very X-ray luminous clusters at z > 0.3, 
thereby more than tripling the number of suc h systems previously 
known. The high-redshift MACS subsample dEbeling et al]|2007D 
comprises 12 clusters a z > 0.5. MACSJ0717.5+3745 is one of 
them. All 12 were observed with the ACIS-I imaging spectrograph 
onboard the Chandra X-ray Observatory. Moderately deep optical 
images covering 30 x 27 arcmin 2 were obtained in five passbands 
(B, V, R, I, z') with the SuprimeCam wide-field imager on the 
Subaru 8.2m Telescope, and supplemented with u-band imaging 
obtained with MegaCam on the Canada France Hawaii Telescope 
(CFHT). Finally, the cores of all clusters in this MACS subsam- 
ple were observed with the Advanced Camera for Surveys (ACS) 
onboard HST in two bands, F555W & F814W, for 4.5ks in both 
bands, as part of programmes GO-09722 and GO- 11560 (PI Ebel- 
ing). 

2.1 HST/ACS Wide-Field Imaging 

A mosaic of images of MACSJ0717.5+3745 and the filamentary 
structure to the South-East was obtained between January 24 and 



Table 2. Overview of groundbased imaging observations of 
MACSJ0717.5+3745. 









Subaru 








CFHT 






B 


V 


R C 


k 


z' 


u* 


J 


K S 


Exposure (hr) 


0.4 


0.6 


0.8 


0.4 


0.5 


1.9 


1.8 


1.7 


Seeing (arcsec) 


0.8 


0.7 


1.0 


0.8 


0.6 


1.0 


0.9 


0.7 



February 9, 2005, with the ACS aboard HST (GO-10420, PI Ebel- 
ing). The 3x6 mosaic consists of images in the F606W and F8 14W 
filters, observed for roughly 2.0 ks and 4.0 ks respectively (1 & 2 
HST orbits). Only 17 of the 18 tiles of the mosaic were covered 
though, since the core of the cluster had been observed already (see 
Tab.Qjfor more details). 

Charge Transfer Inefficiency (CTI), due to radiation damage 
of the ACS CCDs above the Earth's atmosphere, creates spurious 
trails behind objects in HST/ACS images. Since CTI affects galaxy 
photometry, astrometry, and shape measurements, correcting the ef- 
fect is critic al for weak-lensing studies. We apply the algorithm 
proposed by iMassev et al.l j20ld) which operates on the raw data 
and returns individual electrons to the pixels from which they were 
erroneously dragged during readout. Image registration, geometric 
distortion corrections, sky subtraction, cosmic ray rejection, and the 
final combination of the dithered images are then performed usin g 
the standard MULTIDRIZZLE routines jKoekemoer et alj|2002l) . 
MULTIDRIZZLE parameters are set to values optim ised for pre- 
cise galaxy shape measurement dRhodes et al]|2007l) . and output 
images created with a 0.03" pixel grid, compared to the native ACS 
pixel scale of 0.05". 
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Figure 2. Colour-colour diagram (B-R vs R-I) for objects within the 
HST/ACS mosaic of MACSJ0717.5+3745. Grey dots represent all objects 
in the study area. Unlensed galaxies diluting the shear signal are marked by 
different colours: galaxies spectroscopically confirmed as cluster members 
or foreground galaxies (green); galaxies classified as foreground objects be- 
cause of their photometric redshifts (red); and galaxies classified as cluster 
members via photometric redshifts (yellow). The solid black lines delin- 
eate the BRI colour-cut defined for this work to mitigate shear dilution by 
unlensed galaxies. 



2.2 Groundbased Imaging 

MACSJ07 17.5+3745 was observed in the B, V, R c , I c and z' bands 
with t he Suprime-Cam wide -field camera on the Subaru 8.2m tele- 
scope dMivazaki et alj2o"()2l) . These observations are supplemented 
by images in the u* band obtained with the MegaPrime camera on 
the CFHT 3.6m telescope, as well as near-infrared imaging in the 
J and Ks bands obtained with WIRcam on CFHT. Exposure times 
and seeing conditions for these observations are listed in Tab. [2] 
(see also C.-J. Ma, Ph.D. thesis). All data were reduced using stan- 
dard techniques which were, however, adapted to deal with special 
characteri stics of the Supri me-Cam and MegaPrime data; for more 
details see lDonovanl | |2007|) . 

The groundbased imaging data thus obtained are used primar- 
ily to compute photometric redshifts which allow the elimination 
of cluster members and foreground galaxies that would otherwise 
dilute th e shear signal. To this end we use the object catalogue com- 
piled by |MaetalJ<2008l) which we describe briefly in the follow- 
ing. Imaging data from the passbands liste d in Tab. [2] were s eeing- 
matched using the technique described in iKartaltepe etafl <2008l) 
in order to allow a robust estimate of the spectral energy distribu- 
tion (SED) for all objects within the field of view. The object cata- 
logue was then created u sing the SExtractor photometry package 
jBertin & Arnoutsl 19961) in "dual mode" using the R-band image as 
the ref erence detection image. More details are given in lMa et all 
d2008t> . 



Figure 3. Redshift distribution of all galaxies with B, R c , and I c photometry 
from Subaru/SuprimeCam observations that have photometric or spectro- 
scopic redshifts (black histogram). The cyan histogram shows the redshift 
distribution of galaxies classified as background objects using the BRI cri- 
terion illustrated in Fig. [2] 



2.3 Spectroscopic and Photometric Redshifts 

Spectroscopic observations of MACSI07 17.5+3745 (including the 
full length of the filament) were conducted between 2000 and 2008, 
mainly with the DEIMOS spectrograph on the Keck-II 10m tele- 
scope on Mauna Kea, supplemented by observations of the cluster 
core region performed with the LRIS and GMOS spectrographs on 
Keck-I and Gemini-North, respectively. The DEIMOS instrument 
setup combined the 600ZD grating with the GC455 order-blocking 
filter and a central wavelength between 6300 and 7000 A; the expo- 
sure time per MOS (multi-object spectroscopy) mask was typically 
3x1800 s. A total of 18 MOS masks were used in our DEIMOS 
observations; spectra of 1752 unique objects were obtained (65 of 
them with LRIS, and 48 with GMOS), yielding 1079 redshifts, 537 
of them of cluster members. Figure [T] shows the area covered by 
our spectroscopic survey as well as the loci of the targeted galax- 
ies. The data were reduced with the DEIMOS pipeline developed 
by the DEEP2 project. 

Photometric redshifts for galaxies with <2 4.0 were com- 
puted using the adaptive SED- fitting code Le Phare jArnouts et al.l 
1 19991 : lllbert et alj|200d hoO^t) . In addition to employing X 2 opti- 
mization during SED fitting, Le Phare adaptively adjusts the pho- 
tometric zero points by using galaxies with spectroscopic redshifts 
as a training set. This approach reduces the fraction of catastrophic 
errors and also mitigates systematic trends in the differences b e- 
tween spectroscopic and photometric redshifts Gilbert et all 20061) . 

Further details, e.g. concerning the selection of targets for 
spectroscopy or the spectral templates u sed for the dete rmination 
of photometric redshifts, are provided bv lMa et"al] J2008l) . The full 
redshift catalogue as well as an analysis of cluster substructure and 
dynamics as revealed by radial velocities will be presented in Ebel- 
ing et al. (2012, in preparation). 
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Figure 4. As Fig. [2] but for the B-V and u-B colours. The solid black 
line delineates the uBV colour-cut defined for this work to mitigate shear 
dilution by unlensed galaxies. 

3 WEAK GRAVITATIONAL LENSING ANALYSIS 

3.1 The ACS catalogue 

Our weak-lensing analysis is based on shape measurements in the 
ACS/F814W band. Following a method developed for the anal- 
ysis of data obt a ined f or the COSMOS survey and described in 
iLeauthaud et alj J2007f) (hereafter L07) we use the SExtractor 
photometry package dBertin & Arnoutsll 19961) to detect sources in 
our ACS ima ging data in a two-step process. Called the "Hot-Cold" 
technique l lRix et alj2004l L07), it consists of running SExtractor 
twice: first with a configuration optimised for the detection of only 
the brightest objects (the "cold" step), then a second time with a 
configuration optimised for the detection of the faint objects (the 
"hot" step) that contain most of the lensing signal. The resulting 
object catalogue is then cleaned by removing spurious or duplicate 
detections using a semi-automatic algorithm that defines polygonal 
masks around stars or saturated pixels. 

Star-galaxy classification is performed by examining the dis- 
tribution of objects in the magnitude (MAG-AUTO) vs peak 
surface-brightness (MUJV1AX) plane. This diagram allows us to 
separate three classes of objects: galaxies, stars, and any remaining 
spurious detections (i.e., artifacts, hot pixels and residual cosmic 
rays). Finally, the drizzling process introduces pattern-dependent 
correlations between neighbouring pixels which artificially reduces 
the noise level of co-added drizzled images. We apply the remedy 
used by L07 by simply scaling up the no ise level in each pixel b y 
the same constant Fa ~ 0.316, defined bv lCasertano et al.1 12000l) . 

3.2 Foreground, Cluster & Background Galaxy 
Identifications 

Since only galaxies behind the cluster are gravitationally lensed, the 
presence of cluster members and foreground galaxies in our ACS 
catalogue dilutes the observed shear and reduces the significance of 



Figure 5. Redshift distribution of all galaxies with u, B, and V photome- 
try from CFHT/MegaCam and Subaru/SuprimeCam observations that have 
photometric or spectroscopic redshifts (black histogram). The cyan his- 
togram shows the redshift distribution of galaxies classified as background 
objects using the uBV criterion illustrated in Fig. [4] 

all quantities derived from it. Identifying and eliminating as many 
of the contaminating unlensed galaxies is thus critical. 

As a first step, we identify cluster galaxies with the help 
of the ca talogue o f photo metric and spectroscopic redshifts com- 
piled by iMa etalJ d2008f) from groundbased observations of the 
MACSJ07 17.5+3745 field; the limiting magnitude of this cata- 
logue is m Rc = 24. According to lMa & Ebelina ( 1201 Of) , all galaxies 
with spectroscopic redshifts 0.522 < z spec < 0.566 and with photo- 
metric redshifts 0.48 < z p i,ot < 0.61 can be considered to be cluster 
galaxies. An additional criterion can be defined using the photomet- 
ric redshifts derived as described in Sect. 12.31 Taking into account 
the statistical uncertainty of Az = 0.021 of the photometric red- 
shifts, galaxies are defined as cluster members if their photometric 
redshift satisfies the criterion: 

\Zplwl Zch{sler\ ^ (Tphot—zi 

with 

(Tphot-z = (1 + Zdu SKr )kz = 0.036, 

where z p i,ot and Zduster are the photometric redshift of the galaxy 
and the spectroscopic redshift of the galaxy cluster respectively. 
Reflecting the need for a balance between completeness and con- 
tamination, these redshift limits are much more generous than those 
used in conjunction with spectroscopic redshifts, which set the red- 
shift range for cluster membership to 3<x ~ 0.0122. For mor e 
details, see C. - J. Ma (Ph.D. thesis), as well as iMa et al.1 ( 120081) : 
lMa&Ebelingl ( l2010l) . 

In spite of these cuts according to galaxy redshift, the remain- 
ing ACS galaxy sample is most likely still contaminated by fore- 
ground and cluster galaxies, the primary reason being the large 
difference in angular resolution and depth between the ACS and 
Subaru images. The relatively low resolution of the groundbased 
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data causes the Subaru catalogue to be confusion limited and makes 
matching galaxies between the two catalogues difficult, especially 
near the cluster core. As a result, we can assign a redshift to only 
~ 15% of the galaxies in the HST/ACS galaxy catalogue. 

For galaxies without redshifts, we use colour-colour diagrams 
(B-R vs R-I, Fig. |2l and B-V vs u-B, Fig.© to identify fore- 
ground and cluster members. Using galaxies with spectroscopic or 
photometric redshifts from the full photometric Subaru catalogue 
with a magnitude limit of zn R[ . = 25, we identify regions marked 
dominated by unlensed galaxies (foreground galaxies and cluster 
members). In the BRI plane we find B-R < 2.6 (R-I) + 0.05; 
(R-I) > 1.03; or (B-R) < 0.9 to best isolate unlensed galaxies; in 
the UBV plane the most efficient criterion is B-V < -0.5 (u-B) + 
0.85. 

Figures [3] and [5] show the galaxy redshift distributions before 
and after the colour-colour cuts (BRI or uBV) are applied. The re- 
sults of either kind of filtering are similar. The uBV selection is 
more efficient at removing cluster members and foreground galax- 
ies at z < 0.6 (20% remain compared to 30% for the BRI criterion) 
but also erroneously eliminates part of the background galaxy pop- 
ulation. Comparing the convergence maps for both colour-colour 
selection schemes, we find the uBV selection to yield a better de- 
tection of structures in the area surrounding the cluster, indicating 
that suppressing contamination by unlensed galaxies is more im- 
portant than a moderate loss of background galaxies from our final 
catalogue (see Sect.[5]for more details). 

Since the redshift distribution of the background population 
peaks at 0.61 < z < 0.70 (cyan curve in Fig. [5} we assign, in the 
mass modelling phase, a redshift z = 0.65 to background galaxies 
without redshift. 



3.3 Shape measurements of Galaxies & Lensing Cuts 

3.3.1 Theoretical Weak Gravitational Lensing Background 

The shear signal contained in the shapes of lensed background 
galaxies is induced by a given foreground mass distribution. In the 
weak-lensing regime this shear is observed as a statistical deforma- 
tion of background sources. The observed shape of a source galaxy, 
e, is directly related to the lensing-induced shear, y, according to 
the relation : 

£ ^intrinsic ^lensing' 

where e intlinsic is the intrinsic shape of the source galaxy (which 
would be observed in the absence of gravitational lensing), and 



Kid 



y 



1 -K 

Here k is the convergence. In the weak-lensing regime, k <k 1, 
which reduces the relation between the intrinsic and the observed 
shape of a source galaxy to 

£ ^intrinsic 7* 

Assuming galaxies are randomly oriented on the sky, the ellipticity 
of galaxies is an unbiased estimator of the shear, down to a limit 
refer red to as "intrinsic shap e noise", o~ intrinsic (for more details see 
L07; lLeauthaud et alj|201of , hereafter L10). Unavoidable errors in 
the galaxy shape measurement are accounted for by adding them in 
quadrature to the"intrinsic shape noise": 



The shear signal induced on a background source by a given 
foreground mass distribution will depend on the configuration of 



the lens-source system. The convergence k is defined as the dimen- 
sionless surface mass density of the lens: 



(1) 



where 6 is the angular position of the background galaxy, <p is the 
deflection potential, ~L(DolS) is the physical surface mass density 
of the lens, and E cril is the critical surface mass density defined as 



C 2 Dr. 



AtxGDolD, 



Here, D 0L , D os , and D LS represent the angular distances from the 
observer to the lens, from the observer to the source, and from the 
lens to the source, respectively. 

Considering the shear y as a complex number, we define 

r = 7i + '72, 

where y\ = |y| cos 20 and y 2 = M sin 20 are the two components of 
the shear, y, defined previously, and is the orientation angle. With 
this definition, the shear is defined in terms of the derivatives of the 
deflection potential as : 



yi 



-(ip u - yj 22 ), 



with 



<Pij 



72 = </>12 = <P2l, 



Following Kaiser & Squires (1993), the complex shear is re- 
lated to the convergence by: 

k(0) = -- J d 2 G"Re[D{G - 9')y*(6')]. 

Here D(8) is the complex kernel, defined as 

£>(9) 



e]-e\ + ne l e 1 



and Re(x) defines the real part of the complex number x. The aster- 
isk denotes complex conjugation. The last equation shows that the 
surface mass density k{0) of the lens can be reconstructed straight- 
forwardly if the shear y{6) caused by the deflector can be measured 
locally as a function of the angular position 0. 



3.3.2 The RRG method 

To measure the sh ape of galaxies we use the RRG method 
jRhodesetai]|2000h and the pipeline developed by L07. Having 
been developed for the analysis of data obtained from space, the 
RRG method is ideally suited for use with a small, diffraction- 
limited PSF as it decreases the noise on the shear estimators by 
correcting each moment of the PSF linearly, and only dividing them 
at the very end to compute an ellipticity. 

The ACS PSF is not as stable as one might expect from a 
space-based camera. lRhodesetal1 ( l2007l) showed that both the size 
and the ellipticity pattern of the PSF varies considerably on time 
scales of weeks due to telescope 'breathing'. The thermal expan- 
sion and contraction of the telescope alter the distance between the 
primary and the secondary mirrors, inducing a deviation of the ef- 
fective focus and thus from the nominal PSF which becomes larger 
and more elliptical. Using version 6.3 of the Tiny Tim ray-tracing 
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program, iRhodes et all i l2007t) created a grid of simulated PSF im- 
ages at varying focus offsets. By com paring the ell i pticity of ~ 
20 stars in each image to these models, IRhodes et alj (2007) were 
able to determine the effective focus of the images. Tests of this 
algorithm on ACS/WFC images of dense stellar fields confirmed 
that the best-fit effective focus can be repeatedly determined from 
a random sample of 10 stars brighter than m F814W = 23 with an 
rms error of l//m. Once images have been grouped by their effec- 
tive focus position, the few stars in each images can be combined 
into one large catalog. PSF parameters are then inter polated using 
a pol ynomial fit in the usual weak-lensing fashion dMassev et al.1 
2002h. More detail s on the PSF modelling scheme are given in 
Rhodes et af] | |2007|) . 

The RRG method returns three parameters: d, a measure of the 
galaxy size, and, the ellipticity represented by the vector e = (e\ , e2 ) 
defined as follows: 



e 



a 2 -b 2 

a 2 + b 2 
ecos(20) 

e sin(20), 



where a and b are the half-major and half-minor axis of the back- 
ground galaxy, respectively, and if> is the orientation angle of the 
ellipse defined previously. The ellipticity e is then calibrated by a 
factor called shear polarizability, G, to obtain the shear estimator y: 



y 



4 



(2) 



The shear susceptibility factor G is measured from moments of the 
global di stribution of e and other shape parameters of higher or- 
der (see IRhodes et al.|[200oh . The Shear TEsting Program (STEP; 
iMassev et al.l2007l) for COSMOS images showed that G is not con- 
stant but varies as a function of redshift and S/N. To determine G 
for our galaxy sample we use the same definition as the one used 
for the COSMOS weak-lensing catalogue (see L07): 



1.125 + 0.04 arctan 



S/N- 17 



Finally, C, in Eq. [2] is the calibration factor. It was determined 
using a set of simulated images similar to those used by STEP 
faevmans etali 2006; Ma ssev et alJ r2006 ) for COSMOS images, 
and is given by C = (0.86*gQj) (for more details see L07). 



3. 3. 3 Error of the Shear Estimator 

As explained in Sect. 13.3.11 the uncertainty in our shear estimator 
is a combination of intrinsic shape noise and shape measurement 
error: 



i, + 0-; 



measurement' 



where <x? is referred to as shape noise. The shape measurement 
error is determined for each galaxy as a function of size and mag- 
nitude. Applying th e method implemented in the PHOTO pipeline 
dLupton et al.l200ih to analyze data from the Sloan Digital Sky Sur- 
vey, we assume that the optical moments of each object are the 
same as the moments computed for a best-fit Gaussian. Since the 
ellipticity components (which are uncorrelated) are derived from 
the moments, the variances of the ellipticity components can be 
obtained by linearly propagating the covariance matrix of the mo- 
ments. The value of the intrinsic shape noise, <x mtrinsic , is taken to be 
0.27 (for more details see L07, L10). 



In order to optimize the signal-to-noise ratio, we introduce an 
inverse-variance weighting scheme following L10: 

1 

Wy = — . 

°1 

Hence faint small galaxies which have large measurement errors 
are down-weighted with respect to sources that have well measured 
shapes. 



3.3.4 Lensing Cuts 

The last step in constructing the weak-lensing catalogue for the 
MACSJ0717. 5+3745 field consists of applying lensing cuts, i.e., 
to exclude galaxies whose shape parameters are ill-determined and 
will increase the noise in the shear measurement more than they add 
to the shear signal. However, in doing so, we need to take care not 
to introduce any biases. We use three galaxy properties to establish 
the following selection criteria: 

• Their estimated detection significance: 

S _ EEUXJUJTO 45 
N ~ FLUXERRJWTO > ' ' 

where FLUX_AUTO and FLUXERR.AUTO are parameters re- 
turned by S Extractor; 

• Their total ellipticity: 



e = ^e\+e 2 2 < 1; 
• Their size as defined by the RRG d parameter: 
d > 0.13". 

The requirement that the galaxy ellipticity be less than unity 
may appear trivial and superfluous. In practice it is meaningful 
though since the RRG method allows measured ellipticity values 
to be greater than 1 because of noise, although ellipticity is by def- 
inition restricted to e < 1. Because Lenstool prevents ellipticities 
to be larger than 1, we removed the 251 objects with an ellipticity 
greater than unity from the RRG catalogue (2% of the catalogue). 
Serving a similar purpose, the restriction in the RRG size parameter 
d aims to eliminate sources with uncertain shapes. PSF corrections 
become increasingly significant as the size of a galaxy approaches 
that of the PSF, making the intrinsic shape of a galaxy difficult to 
measure. 

Our final weak-lensing catalogue is composed of 10170 back- 
ground galaxies, corresponding to a density of ~ 52 galaxies 
arcmin~ 2 . In addition to applying the aforementioned cuts, and in 
order to ensure an unbiased mass reconstruction in the weak lens- 
ing regime only, we also remove all background galaxies located 
in the multiple-image (strong-lensing) region defined by an ellipse 
aligned with the cluster elongation and with a semi-major axis of 
55", a semi-minor axis of 33". From the resulting mass-map pre- 
sented in Sect. [5] we a-posteriori derived a convergence histogram 
of the pixels outside this region, and found that 90% of them are 
smaller than k = 0.1, with a mean value of k = 0.03. The shear 
and convergence are so weak because the ratio D LS [D s =0.14 
for background galaxies at redshift z me a = 0-65 and the cluster at 
z = 0.54 (see Sect.gJ. 
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4 MASS DISTRIBUTION 

The mass modelling for the entir e MACSJ0717.5+ 3745 field is 
performed using the LENSTOOlQ jjulloetaljfe ocm software, us- 



Jullo& Kneibl [2009) 



ing the adaptive-grid technique developed by . 
and modified by us for weak-lensing mass measurements. Because 
LENSTOOL implements a Bayesian sampler, it provides many 
mass maps fitting the data that can be used to obtain a mean mass 
map and to determine its error. 



4.1 Multi-Scale Grid Method 

We start with the method proposed by ljullo & Kneibl | |2009|) to 
model the cluster mass distribution using gravitational lensing. This 
recipe uses a multi-scale grid of Radial Basis Functions (RBFs) 
with physically motivated profiles and lensing properties. With a 
minimal number of parameters, the grid of RBFs of different sizes 
provides higher resolution and sharper contrast in regions of higher 
density where the number of constraints is generally higher. It is 
well suited to describe irregular mass distributions like the one in- 
vestigated here. 

The initial multi-scale grid is created from a smoothed map 
of the cluster K-band light and is recursively refined in the densest 
regions. In the case of MACSJ0717.5+3745, this method is fully 
adaptive as we want to sample a wide range of masses, from the 
cluster core to the far edge of the HST/ACS field where the fil- 
amentary structure is least dense. Initially, the field of interest is 
limited to a hexagon, centred on the cluster core and split into six 
equilateral triangles (see Fig. 1 in lJullo & Kneibl2009h . This initial 
grid is subsequently refined by applying a splitting criterion that is 
based on the surface density of the light map. Hence, a triangle will 
be split into four smaller triangles if it contains a single pixel that 
exceeds a predefined light-surface-density threshold. 

Once the adaptive grid is set up, RBFs described by Trun- 
cated Isothermal Mass Distributions (TIMD), circular versions 
of truncated Pseud o Isothermal Elliptical Mass Distributions 
(PIEMD) (see e.g. iKassiola & KovneJll993l : iKneib et aljfl996l ; 
iLimousin et al 1 l2005l : lElfasdottir et al I2007F are placed at the grid 
node locations. The analytical expression of the TIMD mass sur- 
face density is given by 



„) 



with 



f(R,r core ,r cat ') ■ 



2G r cu 



+ R 2 



+ R 2 



Hence, / defines the profile, and <xj defines the weight of the 
RBF. This profile is characterized by two changes in slope at ra- 
dius values of r corc and r cut . Within r core , the surface density is ap- 
proximately constant, between r core and r cut , it is isothermal (i.e., 
£ oc r~'), and beyond r cut it falls as X cc r~ 3 . This profile is physi- 
cally motivated and meets the three important criteria of a) featur- 
ing a finite total mass, b) featuring a finite central density, and c) 
being capable of describing extended flat regions, in particular in 
the centre of clusters. 

The RBFs' r core value is set to the size of the smallest nearby 
triangle, and their r cu , parameter is set to 3r core , a scaling that 



1 LENSTOOL is available online: http://lamwws.oamp.fr/lenstool 




Figure 6. Distribution of grid nodes and galaxy-scale potentials, superim- 
posed on the K-band light map of MACSJ0717. 5+3745. Cyan crosses rep- 
resent the location of 468 RBFs at the nodes of the multi-scale grid; magenta 
circles represent galaxy-scale potentials used to describe the contribution of 
individual cluster members. The white dashed line defines the HST/ACS 
field of view. The white cross marks the cluster centre at 07:17:30.025, 
+37:45:18.58 (R.A., Dec). 



ljullo & Kneibl d2009h found to yield an optimal compromise be- 
tween model flexibility and over fitting. We then adapt the technique 
proposed by iDiego et alj J2007h to our multiscale grid model to fit 
the weak-lensing data (k <k 1). Assuming a set of M images and a 
model comprised of N RBFs, the relation between the weights <x 2 ' 
and the 2M components of the shear is given by 











rf 

















A (1,1) 



v(l.M) 



1 (M,N) 



\ (M,N) 



with 



UJ) _ u ui u.j) 
u os 



A™ = 



D> 



« „(>«/) 
"2 ' 



where 



JhJ) 



D 



(d n <I>j(R,j)-d22®i(Rij)), 



and 



In the above, D' LS and D' os are the angular distances from the 
RBF j to the background source i and from the observer to the 
background source i, respectively, and O is the projected gravita- 
tional potential. 7^=1,2 are the two components of the shear, Arf? 2 



is the value of the RBF normalized with cr: 



1, centered on the 



grid node located at 6j, and computed at a radius R = \0 t - 6j\. 
The contribution of this RBF to the predicted shear at location is 
given by the product A^ l2 crl J (see Eq.Q}. More details about this 
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method will be provided in a forthcoming paper (Julio et al. 2012, 
in preparation). 

Figure[6]shows the distribution of the 468 RBFs superimposed 
on the light map used to define the multi-scale grid. The smallest 
RBF has a r core = 26". Note that although Lenstool can perform 
strong lensing, and reduced-shear optimization, here the proposed 
formalism assumes weak-shear with k <k 1 . To get an unbiased es- 
timator, we remove the background galaxies near the cluster core 
from the catalogue, but we keep RBFs in this region as we found it 
improves the mass reconstruction. In order to check our weak-shear 
assumption, we multiply the resulting convergence map by the fac- 
tor 1 - k, and find this minor correction only represents about 5% at 
300 kpc and less than 1% at 500 kpc. It confirms that working with 
a weak-shear assumption in this case is not biasing our results. On 
Fig. [8] the black dashed curve corresponds to the corrected weak- 
shear optimization, while the black curve represents the weak-shear 
optimization profile. Both shows a really good agreement with the 
strong-lensing results (magenta curve). 




Figure 7. Contours of the convergence k in the MACS0717. 5+3745 field 
obtained using the inversion method described by Seitz & Schneider (1995), 
overlaid on the K-band light map. Magenta contours represent the 1 , 2 and 
3tr contours. Cyan crosses mark the position of two galaxy groups (see text 
for details). The white cross marks again the cluster centre (cf. Fig. [6). 



4.2 Cluster Member Galaxies 

Our catalogue o f cluster membe rs is compiled from the photomet- 
ric catalogue of iMa et al.l | |2008|) . Within the HST/ACS study area 
we identify as cluster members 1244 galaxies with spectroscopic 
and/or photometric redshifts within the redshift ranges defined in 
Sect. 13. 21 We use these galaxies' X"s-band luminosities as mass es- 
timators (see below for details). 

All cluster members are included in the lens model in the form 
of truncated PIEMD potentials (see Sect. 14.11 ) with characteristic 
properties scaled according to their K-band luminosity: 



and. 



a- 



r core( L ») • 

,/M l/4 



These scaling relati ons are found to well describe e arly- 
type cluster galaxies (e.g. IWuvts et alj|2004l ; | Fritz et al. 1 120051) un- 
der the assumpti o ns of mass tracing light and the validity of the 
iFaber & Jacksonl J 19761) relation. Since the mass, M, is propor- 
tional to £r?r cut , the above relations ensure M oc L, assuming that 
the mass-to-light ratio is constant for all cluster members. 

To find suitable values of r* ore , r* ut , and a* Q we take advantage 
of the results of the str ong-lensing analysis of MACSJ07 17.5+3745 
recently conducted bv lLimousin e t al. (2012). Their mass model of 
the cluster also included clu ster members, using the scaling rela- 
tions defined above (see also lLimousin et al]|20071) . For a given L* 
luminosity, given by m* = 1 9.16, the mean appare nt magnitude of 
a cluster member in K-band. [Limousin et al .1 J2012t> set all geomet- 
rical galaxy parameters (centre, ellipticity, position angle) to the 
values measured with SExtractor, fixed r* ore at 0.3 kpc, and then 
searched for the values of <Xg and r* ut that yield the best fit. We here 
use the same best-fit values, r* CM = 60 kpc and cr* a = 163 kpc/s, to 
define the potentials of the cluster members. Hence, all parameters 
describing cluster members are fixed in our model; their positions 
are marked by the magenta circles in Fig. [6] 



4.3 Mass Modeling 

The mass reconstruction is conducted using LENSTOOL which 
implements an optimisation method based on a Bayesian M arkov 
Chain Monte Carlo (MCMC) approach d Julio et al J 120071) . We 
chose this approach because we want to propagate as transparently 
as possible errors on the ellipticity into errors on the filament mass 
measurement. This method provides two levels of inference: pa- 
rameter space exploration and model comparison, by means of the 
posterior Probability Density Function (PDF) and the Bayesian ev- 
idence, respectively. 

All of these quantities are related by the Bayes Theorem: 



Pr(0\D, M) oc 



Pr(D|0, M)Pi(0\M) 
Pr(D\M) ' 



where Pr(#|D, M) is the posterior PDF, Pr(D\6, M) is the likelihood 
of a realisation yielding the observed data D given the parameters 
of the model M, and, ~Pi(0\M) is the prior PDF for the parame- 
ters. Pr(D|M) is the probability of obtaining the observed data D 
given the assumed model M, also called the Bayesian evidence. It 
measures the complexity of the model. The posterior PDF will be 
the highest for the set of para meters tha t yield s the best fit and 
is consistent with the prior PDF. Ijullo et all d2007l) implemented an 
annealed Markov Chain to converge progressively from the prior 
PDF to the posterior PDF. 

For the weak-lensing mass mapping, we implement an addi- 
tional level of complexity based on the Gibbs sampling tech nique 
(see Massive Inference in the Bayesys manual. Iskillindl 1 9981) . Ba- 
sically, only the posterior distribution of the most relevant RBFs is 
explored. The number of RBFs to explore is an additional free pa- 
rameter with a Poisson prior. Exploring possible values of this prior 
in simulations, we find that the input mass is well recovered when 
the mean of this prior is set to 2% of the total number of RBFs in 
the model. The weights of the RBFs, a\ , are decomposed into the 
product of a quantum element of weight, q, common to all RBF, 
and a multiplicative factor In order to have positive masses, we 
make £ follow a Poisson prior (case Masslnf=l), and q to follow 
the prior distribution n(q) = q^qex^ 1 ^ 10 , and we fix the initial 
guess q = 10 km 2 /s 2 . The final distribution of cr 2 ' is well approxi- 
mated by a distribution n(q) with qo = 12 km 2 /s 2 , and, we find that 
16.8 ± 4.8 RBFs are necessary on average to reconstruct the mass 
distribution given our data, i.e. about 3.5% of the total number of 
RBFs in the model. This new algorithm is fast, as it can deliver a 
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Radius (h 74 ' kpc) 

Figure 8. Density profiles within the core of MACS J0717. 5+3749. The ma - 
genta curve represents strong-lensing results from iLimousin et al] feOlA ; 
the black curve shows the profile derived from our weak-lensing analysis; 
the black dashed curve shows the weak-lensing profile corrected from the 
weak-shear approximation 

. The grey area marks the region within which multiple images are observed. 



mass map in 20 minutes for abo ut 10,000 galaxies and 468 RBFs 
— the previous algorithm used in |jullo&Kneibl l l2009h was taking 
more than 4 weeks to converge. The new algorithm has been tested 

in parallel on five processors clocked at 2 GHz. 

W e define the likelihood function as (see e.g. lSchneider et all 
ll2000h : 

Pr(D|0)= i-exp 

where D is the vector of ellipticity component values, and 6 is a 
vector of the free parameters <rj. x 1 is the usual goodness-of-fit 
statistic: 



X' 



^^(r,,-2 r ,,(0,)) 2 



(=i j=i 



(3) 



The 2M intrinsic ellipticity components (M is the number of back- 
ground sources), are defined as follows: 



-j = yj(e j )-2y j (0 j ) 



(4) 



and, are assumed to have been drawn from a Gaussian distribution 
with variance defined in Sect. |3.3.T| : 



The normalization factor is given by 



Z, = 



Note the factor of 2 in Eqn.[3]and[4]because of the particular 
definition of the ellipticity in RRG (see Sect. l3.3.2t . The logarith- 




Figure 9. Mass distribution in the core of MACSJ0717. 5+3745. Cyan con- 
tours represent t he ma ss distribution inferred in the strong-lensing study by 
ILimousin et aljjz012l) ; white contours show the mass distribution obtained 
by our weak-lensing analysis; and magenta contours represent the distribu- 
tion of the cluster K-band light. The orange cross marks the cluster center 
adopted in the analysis (cf. Fig.[6). 



mic Bayesian evidence is then given by 



where the average is computed over a set of 10 MCMC realizations 
at any given iteration step /!,, and the integration is performed over 
all iterations Aj from the initial model (A = 0) to the best-fit result 
(A = 1). An increment in dA depends on the variance between the 
10 likelihoods computed at a given iteration, and a convergence 
rate that we set equal to 0. 1 (see Bayesys manual for details). The 
increment gets larger as the algorithm converges towards 1. 



5 RESULTS 

5.1 Kappa Map : Standard mass reconstruction 

We test our catalogue of background galaxies by mapping the con- 
vergence k in the MACSJ0717.5+3745 field. We use the method 
described in Sect. 13.3.11 which is based on the inversion equa- 
tion found by Kai s er & S quires (1993) and developed further by 
ISeitz & Schneider! Jl995h . It shows that the best density recon- 
structions are obtained when the smoothing scale is adapted to the 
strength of the signal. We find a smoothing scale of 3 arcmin to 
provide a good compromise between signal-to-noise considerations 
and map resolution. We estimate the noise directly from the mea- 
sured errors, cr measurement , averaged within the grid cells. 

The resulting K-map is shown in Fig. [7] overlaid on a smoothed 
image of the K-band light from cluster galaxies. We identify two 
substructures south-east of the cluster core whose locations match 
the extent of the filamentary structure seen in the galaxy distribu- 
tion. The first of these, near the cluster core, coincides with the 
beginning of the optical filament; the second falls close to the ap- 
parent end of the filament close to the edge of our study area. A 
simple inversion technique thus already yields a 2 sigma detections 
of parts of the filament. 

The Chandra observation of MACS07 17.5+3745 shows X- 
ray detections of two satellite groups of galaxies embedded in the 
filament (priva te communication H. Ebeling, see also Fig. 1 of 
iMa et alj|2009l , and Fig.[l0]of this paper). Their X -ray coordinates 
(R.A., Dec, J2000) are : i) 07:17:53.618, +37:42:10.46, and, ii) 
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Figure 10. Mass map from our weak-gravitational lensing analysis overlaid on the light map of the mosaic. The two sets of contours show the X-ray surface 
brightness (cyan), and weak-lensing mass (white). The bold white contour corresponds to a density of 1.84 X 10 8 hn M G .kpc~ 2 . The orange cross marks the 
location of the fiducial cluster centre, and the two blue crosses show the positions of the two X-ray detected satellite groups mentioned in Sect. |5. 1 1 The dashed 
cyan lines delineates the edges of the Chandra/ ACIS-I field of view. The yellow box, finally, marks the cluster core region shown in Fig. [9] The magenta line 
emphasises the extent and direction of the large-scale filament. 



07:18:19.074, +37:41:13.14 (cyan crosses in Fig.^. The first of 
these, close to the filament-cluster interface, is detected by the in- 
version technique, but not the second. Conversely, the peak in the 
convergence map in the southeastern corner of our study area is not 
detected in X-rays. 



5.2 Mass Map : Iterative mass reconstruction 

The primary goal of this paper is the detection of the large-scale 
filament with weak-lensing techniques and the characterisation of 
its mass content. To calibrate the mass map of the entire structure 
(i.e., filament + cluster core) and overcome the mass-sheet degen- 
eracy, we compare the weak-lensing mass obtained for the cluster 
core with the tec hnique described in_Sect.[4] with the strong lensing 
(SL) results from lLimousin et aT] d2012l ). 



5.2.1 Cluster Core 

iLimousin et al.l d2012l) inferred a parametric mass model for the 
core of MACSJ0717. 5+3745 using strong-lensing (SL) constraints, 
specifically 15 multiple-image systems identified from multi-color 
data within a single HST/ACS tile. Spectroscopic follow up of 
the lensed features allowed the determination of a well calibrated 
mass model. The cyan contours in Fig. [9] show the mass dis- 
tribution obtained from their analysis. Four mass components, 
associated with the main light components, are needed to sat- 
isfy the observ a tional constraints. The cluster mass reported by 
ILimousin et al. I J2012h for the region covered by a single ACS 
field-of-view is M SL (R < 500 kpc) = (1.06 ± 0.03) 10 15 fCj\ Mq 
where the cluster centre is taken to be at the position quoted be- 
fore (07:17:30.025 +37:45:18.58). This position was adopted as 



the cluster centre because it marks the b arycenter of the Einstein 
ring measured by iMeneghetti et ail ( 12011). Another SL an alysis of 
MACSJ07 17.5+3745 was performed bv lZitrin et alJd2009l) whore- 
port a mass of M SL (R < 350 kpc) = (7.0 ± 0.5) 10 14 h'l M . 

Figure[8]compares the mass density profiles derived from the 
SL analysis (magenta line) and from our WL analysis (black line); 
note the very good agreement. The WL density profile has been cut 
up to ~300 kpc as this region corresponds to the multiple-image 
regime, therefore does not contain any WL constraints. 

In Fig. [8] the SL density profile is extrapolated beyond the 
multiple-image region (the cluster core), while the WL density 
profile is extrapolated into the multiple-image region. The WL 
mass thus obtained for the cluster core is Mwl(S < 500 kpc) = 
(1.04 ± 0.08) IP 15 K,\ M in excel lent agreement with the value 
measured by ILimousin et al. I d2012l) . We also measure M WL (/? < 
350 kpc) = (5.10 ± 0. 54) 10 14 M , in slight disagreement with 
the result reported bv lZitrin et al.l < [2009h . But the excellent agree- 
ment at larger radii between strong- and weak-lensing results with- 
out any adjustments validates our redshift distribution for the back- 
ground sources, and our new reconstruction method. 

Figure [9] compares the mass contours from the two different 
gravitational-lensing analyses for the core of MACSJ0717. 5+3745 
and illustrates the remarkable agreement between these totally in- 
dependent measurements. Our WL analysis requires a bi-modal 
mass concentration which coincide s with the two main co ncentra- 
tions inferred by the SL analysis of ILimousin et al.l d2012l) . 



5.2.2 Detection of a filamentary structure 

Figure [TO] shows the mass map obtained for the whole HST/ACS 
mosaic using the modeling and optimization method described 
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in Sect. f4] Also shown are the X-ray surface brightness as ob- 
served with Chandra/ACIS-I, and the cluster K-band light. A fil- 
amentary structure extending from the south-eastern edge of the 
cluster core to the south-eastern edge of the ACS mosaic area is 
clearly visible in all three datasets. The projected length of the fil- 
ament is measured as 4.5 hz) Mpc, and the mean density in the 
region of the filament is found to be (2.92 ± 0.66) x 10 s h 14 M G 
kpc -2 . The mass concentration detected to the South-East of the 
cluster core corresponds to a low-mass structure comparable to 
a galaxy group and indeed coincides with a poor satellite cluster 
of MACSJ0717.5+3745 that has been independently identified as 
an overdensity of galaxies and as a diffuse peak in the X-ray sur- 
face brightness. The second, more compact, X-ray emitting galaxy 
group (almost due East of the cluster core and outside the filament) 
is also detected. 

To assess the validity of this filament detection, we test the ro- 
bustness of our multi-scale grid optimization method. As a first step 
we create a uniform grid with a RBF of size r C0Te = 21" and 3169 
potentials. We found that using this uniform grid added noise in lo- 
cations where no structures were detected neither in the optical, nor 
in the mass map produced with the multi-scale grid. The Bayesian 
evidence obtained for the uniform grid (pixel size of r con = 21") 
is log(£') = -6589, compared to log(£) = -6586 for the multi- 
scale grid (pixel size of r COTe = 26"). The small difference between 
these two values demonstrates that the reduced number of RBFs in 
our multi-scale grid does not compromise the measurement. This 
is further confirmed by a test using a low-resolution uniform grid 
with a pixel size of r core = 43" comprising 817 RBFs, for which we 
found a Bayesian evidence of log(£) = -6582. All 3 grids equiv- 
alently reproduce the data. Our multi-scale grid model provides an 
intermediate resolution map. 

Having established that the use of an adaptive grid has no 
detrimental effect on the Bayesian evidence, the second test con- 
sisted in increasing the resolution of the multi-scale grid to a pixel 
size of r core = 13", resulting in 2058 potentials. We found that 
the noise increased with the resolution, whereas all real mass con- 
centrations maintained their size. The Bayesian evidence obtained 
with this high-resolution multi-scale grid is log(£') = -6608. A 
WL mass reconstruction with the high-resolution multi-scale grid 
recovers the density profile obtained with the multi-scale grid us- 
ing r mre = 26" pixels. However, in view of the large difference 
between these two evidence values, it is clear that such an increase 
in resolution is not required by the data. 



6 PROPERTIES OF THE FILAMENT 

In this Section, we first summarise briefly the results of previous 
studies of the MACSJ0717.5+3745 filament, and then discuss the 
filament's properties as derived from our weak-lensing analysis. 

The detection of a filamentary structure i n the field of the clus - 
ter MACSJ0717.5+3745 was first reported bv lEbeling etafl J2004h . 
who discovered a pronounced overdensity of galaxies with V-R 
colours close to the cluster red sequence, extending over ~ 4/j^J 
Mpc at the cluster redshift. Extensive spectroscopic follow-up of 
over 300 of these galaxies in a region covering both the cluster and 
the filament confirmed that the entire structure is located at the clus- 
ter redshift, z = 0.545. This direct detection of an extended large- 
scale filament connected to a massive, distant cluster lent strong 
support to predictions from theoretical models and numerical sim- 
ulations of structure formation in a hierarchical scenario, according 
to which large-scale filaments funnel matter toward galaxy clusters 
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Figure 11. Profile of the mass surface density of the cluster/filament com- 
plex as a function of distance from the cluster core. The three curves show 
the mass surface density measured in our weak-lensing analysis across the 
entire field (black), for the cluster without the filament (cyan), and for the 
inner part of the filament (orange). The two green lines represent two dif- 
ferent slopes: R~ l and R~ 2 . See text for details. 

inhabiting the nodes of the cosmic web. Ta king advantage o f the 
opportunity provided by this special target, iMa et a D d2008h pur- 
sued a wide-field spectroscopic analysis of the galaxy population 
of both the cluster and the filament. Along the filament, they found 
a significant offset in the average redshift of galaxies, correspond- 
ing to ~ 630 km s~' in velocity, as well as a decrease in velocity 
dispersion from the c ore of the cluster to the end of the filament. 
IMa & Ebelingl d2010h studied the morphology of cluster members 
across a 100 fc 2 Mpc 2 area to investigate the effect of cluster envi- 
ronment on star formation using HST/ACS data. They explored the 
relation between galaxy color and density. They showed that the 
1-D density profile of MACSJ07 17.5+3745 cluster members fol- 
lowed the filament from the core to beyond the virial radius falling, 
then slightly rising and flattening at ~ 2/i^ 4 ' Mpc. 

To physically characterise the large-scale filament, we 1) anal- 
yse its density profile (Fig.ll It; 2) study the variation of its diame- 
ter as a function of distance to the cluster core (Fig,ll2t; 3) derive a 
schematic picture of its geometry and orientation along the line of 
sight (Figs. IT31&I 14t: and 4) attempt to measure the mass in stars 
within the entire structure (Fig.|15ll. 

6.1 Density Profiles 

Figure QT] shows the density profile across the entire study area as 
a function of distance from the cluster centre as defined and shown 
in Figs. [6] [7] [9] and[T0] Also shown are the profiles of the filament 
only (obtained by only considering the WL mass along the magenta 
line and within the bold white contour in Fig.llOt. and of the cluster. 
The cluster density profile is cut at 2 ft^'Mpc from the cluster core 
as further, the ACS field is dominated by the filament. 

As is apparent from both Figs. 1 101 and 1111 the filament con- 
nects to the cluster at a distance of about i.Sh^l Mpc from the clus- 
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Figure 12. Diameter of the filamentary structure as a function of distance 
from the cluster centre. The regime of the cluster proper and the segment 
containing the embedded X-ray bright group of galaxies (blue cross in 
Fig. 1101 are marked. The density threshold used to define the filamentary 
structure is 1.84 X 10 s /174 M G .kpc~ 2 . 

ter core, causing the density profile to flatten from this distance 
outward. The mean density in this region is (2.92 ± 0.66) x 10 s /174 
M kpc~ 2 . The embedded galaxy group detected also in the clus- 
ter light and at X-ray wavelengths is responsible for the peak in the 
density profile at a distance of about 2hjl Mpc from the cluster cen- 
tre. The dip in the filament profile at ~2.8 Mpc from the cluster 
centre coincides with a dip visible also in projection in Fig,|10|(see 
Sect. 16.31 for more details on the 3D geometry of the large-scale 
filament). The filament clearly dominates the radial density pro- 
file from ~2 hj^ Mpc from the cluster centre until the edge of the 
HST/ACS field. 

Fitting any of these profiles with parametric models (NFW, 
SIS, etc) is not meaningful in view of the complexity and obvious 
non-sphericity of these structures. We note, however, that the over- 
all density profile decreases as R 2 until about 2 h^l Mpc, where the 
beginning of the filament causes a dramatic flattening of the den- 
sity profile. For illustrative purposes, Fig. [TT] shows two different 
slopes, /J -1 and R 2 . 

6.2 Filament Size 

We define the region of the filament on the mass map by a den- 
sity threshold of 1.84 x 10 s h 14 M G .kpc~ 2 which corresponds to the 
3 sigma detection contour of the lensing mass reconstruction, and 
then define the central axis of the filament with respect to this con- 
tour (magenta line in Fig. |10|. Under the simplifying assumption 
that the filament cross-section is spherical (i.e., the filament resem- 
bles a cylinder of variable radius), the filament diameter is given by 
the perpendicular distance of the contour to this axis. The uncer- 
tainty of this measurement is assessed by repeated this procedure 
for every mass map produced by LENSTOOL, and adopting the 
standard deviation of the results as the error of the filament diame- 
ter. 
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Figure 13. Schematic presentation of the geometry of the large-scale fila- 
ment along the line of sight, as deduced from radial- velocity measurements. 
Over the first 4 Mpc (in projection) the filament's average inclination angle 
is found to be a ~ 75. Beyond this projected distance, the filament curves 
increasingly toward the plane of the sky, eventually rendering its projected 
mass surface density too low to be detectable. In this sketch, radial veloci- 
ties as observed are split into two components: peculiar velocities and Hub- 
ble flow, all expressed relative to the radial velocity of the cluster core. 
(Adapted from Ebeling et al, 2012.) 

Figure [T2l shows the results thus obtained as a function of dis- 
tance from the cluster centre. We find the filament's diameter to de- 
crease with increasing distance from the cluster centre, albeit with 
mild local variations. Recognizable features are the local peak at 
1.85 fijl Mpc from the cluster centre where the presence of an 
embedded group of galaxies (marked by a blue cross in Fig. ITOt 
combined with a second, X-ray faint mass concentration causes the 
filament to widen to (2.79 ± 0.45) Mpc. The following narrow- 
ing of the filament at a cluster-centric distance of about 2.8 Kj\ Mpc 
coincides with a dip in the filament density (cf. Sect. 16. U . Two ad- 
ditional increases in size (neither of them detected in X-rays) are 
observed at larger distances of approximately 3.6 and 4.7 Mpc 
and can again be matched to local peaks in the filament density and 
projected WL mass density (Figs.QT]and[T0] respectively). 



6.3 3D-properties of the large-scale filament 

In order to derive intrinsic properties of the filament from the 
observables measured in projection, we need to know the three- 
dimensional orientation and geometry of the filament and, specif- 
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Figure 14. Mass density of the filament in units of p cr i t as a function of the 
orientation angle a, The red diamond highlights the density for the adopted 
average inclination angle of 75 deg. 



Figure 15. Density profile of the stellar component in the whole HST/ACS 
field of MACSJ07 17.5+3745. The stellar mass is obtained using the 
Arnouts et al. (2007) relation for a red galaxy population, and an Salpeter 
IMF. The mean M./Lj- ratio is equal to 0.94 ± 0.41 for our galaxy popula- 



tion. The two green curves represent two different slopes: R 1 & R 
and dark green respectively. 



, light 



ically, its inclination with respect to the plane of the sky. Ebeling 
et al. (2012, in preparation) derive a model of the entire complex 
by comparing the measured variation in the mean radial velocity 
of galaxies along the filament axis to expectations from Hubble- 
flow velocities, as well as with the predictions of p eculiar veloc- 
ities within filaments from numerical simulati ons dColberg et al .1 
l2005l : ICuesta et alj|2008l : ICeccarelli et alj|201ll) . A self-consistent 
description is obtained for an average inclination angle of the fil- 
ament with respect to the plane of the sky of about 75 degrees 
(Fig.ll3t. At this steep inclination, the three-dimensional length of 
the filament reaches 18 h -1 Mpc within our study region. 

Using the weak-lensing mass map and the above model of the 
3D geometry and orientation of the filament, we can constrain the 
filament's intrinsic density. Adopting again a cylindrical geometry 
and the average projected values from Figs.[TT]and[T2]we obtain a 
mass density of payment = (3-13 ± 0.71) x 10 13 h 2 4 M G Mpc -3 , or 

Pfllamcnt = (206±46)p cri , 

in units of the critical density of the Universe. 

We stress that this is an average value: a more complex density 
distribution along the filament is not only possibly but also physi- 
cally plausible. The density of the filament depends strongly on its 
inclination angle a which is reasonably well constrained on average 
but is uncertain at large cluster-centric distances (see Fig. 1 13l>. Fig- 
ure [14] shows how sensitively the deduced density of the filament 
depends on a. In addition, embedded mass concentrations like the 
X-ray bright group of galaxies highlighted in Fig.[K)]and clearly de- 
tected also in the mass surface density profile (Fig. II lb will cause 
local variations in the filament's density. A more comprehensive 
discussion of dynamical considerations regarding the 3D geome- 
try of the cluster-filament complex is provided in Ebeling et al. (in 
preparation). 



6.4 Stellar Mass Fraction 

To compute the stellar mass, M„, across our study region we use 
the re lation log(M,/L K ) = az + b established by lArnouts et al .1 



2007) for quiescent (red) galaxies in the the VVDS sample 



Le Fevre et al .1 120051) and adopting a Salpeter initial mass func- 
tion (IMF). Here L K is the galaxy's luminosity in the K band, z is 
its redshift, and the parameters a and b are given by 

a = -0.18 ±0.03 b = -0.05 ± 0.03. 

We apply this relation to the K-band magnitudes of all cluster mem- 
bers, defined using the redshift criteria presented in Sect.[3] The K- 
band magnitude limit of our catalogue is 23.1. The resulting pro- 
jected mass density in stars is shown in Fig. [15] as a function of 
distance from the cluster centre. It decreases in proportion to the to- 
tal projected mass density depicted in Fig.QT] as illustrated by the 
two different slopes, S ~ R~ l and S ~ R 2 , shown in both figures. 
Hence the measured mass-to-light ratio is approximately constant 
across the study area, with a mean value of {M,/L K ) = 0.94±0.41. 

T o compare our result with those obtained by Leauthaud et al.l 
j2012l) for COSMOS data we need to adjust our measurement to 
account for the different IMF used by these authors. Applying a 
shift of 0.25 dex to our masses to convert from a Salpeter IMF to a 
Chabrier IMF, we find <M,/L K ktabriei- = 0- 73 ± 0.22 for qu i escen t 
galaxies at z ~ 0.5, in good agreement with lLeauthaud et all d20 1 2|) 



The fraction of the total mass in stars, /„, i.e., the ratio be- 
tween the stellar mass and the total mass of the cluster derived 
from our weak-lensing analysis, is /. = (1.3 ± 0.4)% and /„ = 
(0.9±0.2)% for a Salpeter and a Chabrier IMF, respectively. The lat- 
ter value is sligh t ly low er than that obtained for COSMOS data by 
iLeauthaud et al.l | |2012|) . possibly because of different limiting K- 
band magnitudes or differences in the galaxy environments probed 
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(the COSMOS study was conducted for groups with a halo mass 
comprised between 10" and 10 14 M Q and extrapolated to halos of 
~ 10 15 /r 1 Mo). 

Finally, we compute the total stellar mass within our study area and 
findM, = (8.62±0.24)xl0 13 M G and M, = (6.71±0.19)x 10 13 M G 
for a Salpeter and a Chabrier IMF, respectively. 



7 SUMMARY AND DISCUSSION 

We present the results of a weak-gravitational lensing analysis of 
the massive galaxy cluster MACSJ0717. 5+3745 and its large-scale 
filament, based on a mosaic composed of 18 HST/ACS images 
covering an area of approximately 10 x 20 arcmin 2 . Our mass re- 
const ruction method uses RRG shape measurements l iRhodes et al.l 
120071) ; a multi-scale adaptive grid designed to follow the struc- 
tures' K-band light and including galaxy-size potentials to account 
for cluster members; and the LENSTOOL software package, im- 
proved by the implementation of a Bayesian MCMC optimisation 
method that allows the propagation of measurement uncertainties 
into errors on the filament mass. As a critical step of the analysis, 
we use spectroscopic and photometric redshifts, as well as colour- 
colour cuts, all based on groundbased observations, to eliminate 
foreground galaxies and cluster members, thereby reducing dilu- 
tion of the shear signal from unlensed galaxies. 

A simple converge nce map of the study are a, obtained with 
the inversion method of ISeitz & SchneideJ Jl995l) , already allows 
the detection of the cluster core (at more than 6cr significance) and 
of two extended mass concentrations (at 2 to 3<x significance) near 
the beginning and (apparent) end of the filament. 

The fully optimised weak-lensing mass model yields the sur- 
face mass density shown in Fig. [TO] Its validity is confirmed by the 
excellent agreement between the mass of the cluster core measured 
by us, M WL (R < 500 kpc) = (1.04 ± 0.08) x 10 15 tq\ Mp, and th e 
one obtained in a strong-lensing analysis bv lLimousin et al.l d2012h . 
M SL (R < 500 kpc) = (1.06 ± 0.03) x 10 15 Kj\ M Q . 

Based on our weak-lensing mass reconstruction, we report the 
first unambiguous detection of a large-scale filament fueling the 
growth of a massive galaxy cluster at a node of the Cosmic Web. 
The projected length of the filament is approximately 4.5 tcl Mpc, 
and its mean mass surface density (2.92 ± 0.66) x 10 8 hj4 M Q kpc -2 . 
We measure the width and mass surface density of the filament as 
a function of distance from the cluster centre, and find both to de- 
crease, albeit with local variations due to at least four embedded 
mild mass concentrations. One of these, at a projected distance of 
1.85 hjl Mpc from the cluster centre, coincides with an X-ray de- 
tected galaxy group. The filament is found to narrow at a cluster- 
centric distance of approximately 3.6 Kj\ Mpc as it curves south 
and, most likely, also recedes from us. We find the cluster's mass 
surface density to decrease as r~ 2 until the onset of the filament 
flattens the profile dramatically. 

Following the analysis by Ebeling et al. (2012, in preparation) 
we adopt an average inclination angle with respect to the plane of 
the sky of 75° for the majority of the filament's length. For this sce- 
nario and under the simplifying assumption of a cylindrical cross- 
section we obtain estimates of the filament's three-dimensional 
length and average mass density in units of the Universe's criti- 
cal density of 1 8/774 Mpc and (206 ± 46)p cr jt, respectively. How- 
ever, additional systematic uncertainties enter since either quantity 
is sensitive to the adopted inclination angle. These values lie at the 
high end of the range predicted from numerical simulations (e.g., 



IColberg et al]|2005h . which is not unexpected given the extreme 
mass of MACSJ07 17.5+3745. 

The lensing surface mass density and the spectroscopic red- 
shift distribution suggest the consistent picture of an elongated 
structure at the redshift of the cluster. The galaxy distribution along 
the filam ent appears to be h omogeneous, and at the cluster red- 
shift (see lEbeling et alj|2004h . This motivates our conclusion of an 
unambiguous detection of a large scale filament, and not a super- 
position of galaxy groups projected on the plane of the sky. 

Finally, we measure the stellar mass fraction in the entire 
MACSJ07 17.5+3745 field, using as a proxy the K-band luminosity 
of galaxies with redshifts consistent with that of the cluster-filament 
complex. We find /. = (1.3 ± 0.4)% and /» = (0.9 ± 0.2)% for a 
Salpeter and Chabrier IMF, respectively, in good agreement with 
previ ous results in the fields of massive clusters dLeaufhaud et al .1 
120121) . 

Our results show that, if shear dilution by unlensed galax- 
ies can be efficiently suppressed, weak-lensing studies of mas- 
sive clusters are capable of detecting and mapping the complex 
mass distribution at the vertices of th e cosmic web . Confirming 
results of numerical sim ulations (e.g. IColbergetal]|l999ll2005t 
lOguri & Hamanal 1201 lh , our weak-lensing mass reconstruction 
shows that the contribution from large-scale filaments can be sig- 
nificant and needs to be taken into account in the modelling of mass 
density profiles. Expanding this kind of investigation to HST/ACS- 
based weak-lensing studies of other MACS clusters will allow us to 
constrain the properties of large-scale filaments and the dynamics 
of cluster growth on a sound statistical basis. 
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